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Abstract 

We discuss a simple phenomenological Landau theory of phase transitions with two coupled 
single-component order parameters and compare the results with available experimental data. 
The model corresponds to the case of a ferroic system, in which ferromagnetic and ferroelectric 
transitions originally occur at temperatures Tm and TV, respectively. For Ty > Tm the magne- 
toelectric coupling strongly renormalizes the magnetic transition temperature, Tm — > Trm (with 
Trm » Tm), as well as generates an additional anomaly in ferroelectric subsystem Trm- Full 
susceptibility tensor has also been determined. The concept of Arrot plot is replaced by the Arrot 
planes which appear when both types of order coexist. The results are in good overall agreement 
with experimental data for the ferroelectromagnetic BiMnOs. We also estimate the contribution of 
Gaussian fluctuations of both order parameters, that lead to corrections to the mean- field specific 
heat. Those corrections are still insufficient even though other quantities agree quite well with 
experiment. We calculate the temperature dependence of the coherence length for both types of 
order as well. 
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I. INTRODUCTION 



Multiferroics are materials, in which at least two types of order coexist. For example, 
(anti)ferroelectricity and (anti)ferromagnetism can take place simultaneously [HE]- Multi- 
ferroicity of frustrated magnets, in which magnetism and ferroelectricity coexist with gigantic 
magnetoelectric coupling, has attracted an interest due to challenges to many-body theory, 
as well as by discoveries of new phenomena with a promice for potential applications as 
transducers, actuators, and sensors [3H5]. Quite a few of these multiferroics are manganites, 
in which the magnitude of the spin of the Mn+ 3 ion is large and hence may be treated 
semiclassically. 

The perovskites like AMnC>3 has been widely studied in this context due to the coexistence 
of ferroelectric and magnetic order in some of them [6HSj- In BiMn0 3 the nature of the A 3+ 
ion is central to determine the structural, ferroelectric, and magnetic properties of this 
system [9] . With the help of first principle calculations, Hill et al. [10] provided the reasons 
why we observe so few ferroelectric magnets and predicted the existence of ferroelectricity in 
BiMnCV In ferroelectrics such as SrTiC>3 it is usually driven by a hybridization of empty 3d 
transition metal orbitals with occupied 2p orbitals of the octahedrally coordinated oxygen 
ions. The appearance of magnetic moment in turn, requires partial occupancy of the 3c? 
orbitals. So, in materials such as BiMnC>3, the coexistence of Mn 3+ ions (3d 4 configuration) 
with 6s 2 lone electron pairs due to the Bi 3+ ions can lead to the coexistence of magnetic 
order with electric polarization at low temperature [TT] . 

BiMnC>3 has been the subject of considerable interest, mainly due to its structural sim- 
plicity. From magnetic point of view, Mn 3+ ion in this case has magnetic moment of 3.6/is, 
a value close to ground-state induced by the Hund's rule coupling. In the octahedral 
environment the electronic configuration is t\ g e g . From electrical point of view BiMnC>3 is 
an insulator [T2] . 

A detailed structural study of BiMn0 3 suggests [T31 E] that the material has a highly 
distorted perovskite structure (centrosymmetric space group C2/c) which is incompatible 
with the existence of ferroelectricity. In spite of this, there is an experimental evidence of 
ferroelectric order in a BiMnC>3 polycristaline samples [131 HH1 US] . The discussion concerning 
the origin of the ferroelectricity in this compound is still controversial. 

We present a simple phenomenological approach useful in describing systems like BiMnC>3 
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with two coupled order parameters, as well as compare the results with experiment. This is to 
show to what extent a simple Landau-type approach can account for the experimental results 
in a quantitative manner. Our task is related to an even more basic question to what extent 
ferroelectric order (appearing first at much higher temperature Tf) suppresses the magnetic 
fluctuations near the corresponding transition temperature Tu << Tf. In such situation the 
mean- field description of the magnetic phase should be at least semiquantitatively correct, 
as we demonstrate below. 

The structure of the paper is as follows. In Sections II and III we formulate the Landau- 
theory with two single-component, spatially homogeneous order parameters P and M (rep- 
resenting the ferroelectric and the ferromagnetic types of order, respectively) coupled via 
a phenomenological term of the type — I7I (PM) 2 . This coupling leads to a renormalized 
magnetic transition temperature. In Section III we also compare the results obtained with 
the data available for BiMn03, as well as introduce a new concept of Arrot plates, which 
can be used in systems with two coupled order parameters in a ordered state. In Section 
IV we include Gaussian fluctuations for the coupled system, and subsequently compare the 
results with experimental data concerning the temperature dependence of the specific heat 
near the low-temperature magnetic transition. We summarize our results and provide an 
outlook in Section V. 



II. MEAN FIELD APPROXIMATION: LANDAU APPROACH 

The way to describe the coupling between magnetism and dielectricity in multiferroics 
was proposed by Smolenskii [17] . who explained the origin of the anomaly in the dielectric 
constant in a ferroelectromagnet within the framework of Landau theory of second-order 
phase transitions. In the simplest case, the Landau free energy for a system with two 
coupled order parameters, P and M can be written as: 

F(P, M, T) = F (T) + a o( T - T M) M 2 + h -M A (1) 

ao(T-T f ) 2 4 
2 4 
+ |(PM) 2 -PE a -MH a , 

where E a and H a are the applied electric and magnetic fields, respectively. Recently it 
was shown that this kind of model can discribe the phase transitions on a scale-free network 
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as well [18] . As one can see, the system is characterized by two bare transition temperatures, 
Tf and Tm representing ferroelectric and ferromagnetic transitions, respectively. The form 
of the magnetoelectic coupling term oc 'y(PM) 2 can be explained by using general symmetry 
arguments. Namely, the onset of ferroelectric order requires the breaking of spatial inversion 
symmetry, whereas the appearance of a spontaneous magnetization is connected with the 
breakdown of time reversal symmetry. The coupling term in ^ obeys those two conditions 
and allows for ferroic order with simultaneous nonzero P and M. 



In further calculations we use a dimensionless form of the free energy expansion (|B5|) 
which is obtained from (2) by dividing both sides of (2) by the constant value 
effect, we obtain: 



n 2 T 2 



AF 









if 

V-'M 






+ 








1 


+ 


-7m 



I I M 2 + -M 4 
4 



Tf_ 



T 



l\P 2 + -l 



Tm 



p4 



Tm 



Pe - Mh, 



where: AF 



o 2 T 2 
"0 1 M 



^[F(P,M,T)-F (T)], h 
m = 



i 



l ± M 

M 
m ' 



Ha = ShH a , I 



b ' 



Pa 



Po 



a T, 



f 



7n 



/3a ' 



a 2 b 



(2) 



Within this formulation 7jr-,l , and 7 m are the material parameters. In this case, the bare 
magnetic transition temperature Tm is assumed as always nonzero, whereas the ferroelectric 
temperature is Tf > 0. 



Taking derivatives of (B5) with respect to M and P, respectively, we obtain the following 



system of nonlinear equations: 
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(T/T M - 1)M + M 3 + lm ^-MP 2 -h = 0, 
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(3) 
(4) 



In the limit of zero electric field, we obtain the following relation between P and M from 

P = ±Jl-^-^M 2 . (5) 
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Substituting this expression to (3) we have the equation for magnetization in the usual form: 

A(T)M + BM 3 -h = 0, (6) 

with 

A = A(T) = - 7m ) + lm lL - 1, (7) 

J-M J-M 

and 



B = 1 



/ rr. 

T 



(8) 



One sees that the effective magnetic transition temperature Trm is renormalized by the 
magnetoelectrical coupling. Explicitly, since the renormalized transition temperature is 
determined from the condition A(Trm) — 0, this yields: 



1 - 7r, 
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(9) 



This is one of the interesting results. Namely, the renormalization is strong because of 
negative value of coupling constant 7 m . Furthermore, the renormalization of Tf does not 
appear if Tf > T M . In the case of Tf = T M the coupling would not change the critical 
temperatures and no renormalization would occur. 

The explicit stable solutions of equation uQu for M in the case e = are: 
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(10) 



(v / 3v / 4A 3 B 3 +27B 4 /i 2 -9B 2 /i) 

We apply the solution obtained above to the discussion of selected magnetic and dielectric 
properties of BiMnC>3. 



III. APPLICATION TO BiMn0 3 
A. Magnetic properties 

To visualize the influence of the magnetoelectric coupling on the magnetic properties of 
BiMnC>3, we fitted the temperature and applied magnetic field dependences of the magne- 
tization, based on the data of Kimura [16] and Chiba [19] . In Fig. [l] we display the fitted 
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M(H a ;T) curves near the critical temperature (Trm ~ 100K). In the inset we plot the 
values of A(T) obtained from the fitting: it is indeed a linear function of T, as obtained 
in ([7]). By taking the value Tf = 760K [T6] El we have obtained the renormahzed value 
of Trm ~ 100. 5K, the bare Curie temperature Tm = 2.28K, the magnetoelectric coupling 
constant 7 m ~ —0.15, and I ~ 0.024. One sees that the renormalization Trm/Tm is very 
large with the increase caused by the negative sign of the coupling constant 7. We consider 
the coupling in case of BiMn03 to be large because |7 m | ~ \fl and as it can be seen from 



Eqs. (19-20) it directly effects the rapid increase of electric polarization and magnetization 
in ordered state and is the greatest possible value for the coupling. 

In Fig. [2] we fitted the temperature dependence of magnetization [19] in two ways: first 
(dotted line), by taking the averaged values from Table Jl] and second, by a direct fitting i.e. 
changing slightly the averaged values (dashed line). Those slight changes are justified, as 
they are within statistical error. One should mention that data used in Fig. [I] was taken for 
a different sample to that of Fig. [2j 

From these two figures one sees, that the overall behavior of the magnetization near Trm 
is well reproduced by the mean field approach, particularly for T — > Trm and above. Hence, 
we parametrize the dielectric and magnetoelectric susceptibility components in the same 
manner next. 
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FIG. 1: (Color online). Isothermal magnetization as function of a magnetic field at various tem- 
peratures for BiMnOa [16] (solid line) and that from Eq. ^ (dashed lines). Inset: A(T) values 
for the temperatures marked. The fitting parameters are listed in Table]!] 
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FIG. 2: (Color online). Temperature variation of magnetization of BiMnOs measured at IT [p5] 
(solid line) and the fitted solution of ([6]). Dotted line: The averaged (cf. Table Jl]) fitting parameters 
with A(T) = -50.593 + 0.503T, B = 0.062, m = 0.161 /ie/Mn site, and s h = 26.5 T~ x . A better 
fit (dashed line) can be obtained by a slightly different set of parameters: A(T) = —49. 926+0. 507T, 
B = 0.059, m = 0.208 /i B /Mn site and s h = 24.6 T~ l . 

TABLE I: Fitting parameters obtained from the results of Fig. 1. 



T[K] 


A 


B 


mo [Ms/Mn site ] 


Sh [1/T] 


130 


15.079(8) 


0.0623(4) 


0.1601(8) 


26.509(2) 


110 


4.12(63) 


0.06(233) 


0.16(019) 


26.5(092) 


100 


0.00061 


0.06(047) 


0.16(513) 


26.5(346) 


90 


-5.1(848) 


0.06(234) 


0.16(018) 


26.5(093) 



B. Susceptibility tensor 

Kimura et al. reported [TB] that with increasing temperature the magnitude of the 
isothermal magnetocapacitance increases and exibits a maximum around Trm- Whereas 
upon further increase (above Trm) the magnetocapacitance subsequently decreases. The 
authors claim that this phenomenon arises from the magnetization rotation in magnetic 
domains. Such behavior can also be obtained from the simple Landau approach introduced 
here without involving any domain formation. On application of external fields the system 
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response in the ferromagnetoelectric state (T < Trm) is described by the tensor: 

X = I _ . I ' ( U ) 

Xme Xm 

with: 

<9m _ _ ap _ _ <9m _ _ ap _ _ 

Xmi ~^g~ Xme; Xem- \ / 

By assuming e 7^ 0, we obtain the following equations for M and P: 
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After differentiating (13) and (14) with respect to both e and h, we obtain a system of 
linear equations for the susceptibility components in the form: 



^Xe + C'Xme = 1, 
^Xem + C'Xm = 0, 
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Therefore the solution of (15) takes the form 
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In the zero field case and for temperature T < Trmi the corresponding expressions for 
magnetization and polarization are: 
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FIG. 3: Polarization as a function of temperature for selected values of applied magnetic field. The 
parameters are the same as those used previously. 

In Fig. [3] we plot the temperature dependence of the polarization in the vicinity of 
Trm, i-e. at temperatures T < Tf. One sees that even though the ferroelectric ordering 
temperature is well above Trm, the weaker-scale magnetic interaction significantly enhances 
the polarization. This enhancement is also present on application of a magnetic field. This 
is a second (in addition to the renormalization of Tm) important cross-effect correlating 
magnetic and electric properties in these systems. One should also note that the electric 
polarization increases with increasing magnetic field, as one may expect from the negative 
sign of the magnetoelectric coupling. Unfortunately, no experimental results are available 
to us to confront our findings with those for BiMnOs. 
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We now turn to the analysis of the susceptibility-tensor components. Substituting the 



values obtained above for the magnetization and polarization into (17), we obtain an explicit 
expression for the zero field electric susceptibility, namely 



£ e(0 ) = | l *T f {l(T f -T)+ lm( T-T M) ) f° T T < TrM > 



After calculating x e (h), using the general solutions of (10), we plot: 
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FIG. 4: (Color online). Field-induced change in the dielectric constant as a function of an applied 
magnetic field for selected temperatures. We assume that ((3alT M )/(4ira bT f ) = 0.0008. 

shown in Fig. |4j The curves obtained reflect the corresponding experimental data well 
[IB"] . Though one should note, that the actual change of the dielectric constant is not as 
rapid as in our mean field approach. Also the rapid trend upward of the curves appears 
here above 100K, whereas the respective changes of the data appears only above 110K [16J. 
This difference is a clear sign of the nonzero value of the magnetization due to short-range 
correlation. This type of crossover behavior above Trm will appear in the specific heat data, 
as discussed in the next Section. 
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Experimentaly, the susceptibility is a linear function of the squared magnetization for 
BiMnC>3 [TB] . In some papers |16| |2"U] this is rationalized on the basis of Landau-theory yet 
this it is not the case. In fact, in this phenomenological approach the inverse susceptibility is 
a linear function of the squared magnetization and it comes about from the renormalization 
of dielectric constant by the coupling. Explicitly this may be written [21] AF « (cxq(T — 
Tf) /2 + 7/2 M 2 )P 2 + ... = Xe 1 P 2 + ■■■ ■ I n Fig- [5] we plot the predicted Landau theory value 
of Ae/e(0) as a function of the squared magnetization. The dependence is none linear. In 
the inset of Fig. [5] we show the inverse dielectric susceptibility as a function of the squared 
magnetization to be linear within the framework of Landau theory. 
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FIG. 5: Field induced change in relative dielectric constant as a function of the square of the 
magnetization at 100 K. We assume that (0a§Tjif)/(4.7rao&T/) = 0.0008. Inset: Depend ence of 
the inverse dielectric susceptibility vs. M 2 . The data can be parametrized by the straight line: 
Xe(0)/xe(h) = 1,233 + 1,071M 2 . 



In case of Xme = we calculate the inverse susceptibilities from the Landau functional as 
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in which P and M may be evaluated from Eqs. ^ and (10). In Fig. [fj] and [7] we plot the 
calculated temperature dependences of the susceptibilities mentioned above. The magne- 
toelectric coupling causes the enhancement of polarization around Trm and a suppression 
of the dielectric constant cf. Fig. |6j One can see that at zero applied magnetic field % e 
decreases stepwise at Trm whereas it is gradually suppressed with increasing field. Xe(T) 
exhibits a trend observed experimentally with increasing H a , but the calculated changes 
are too large. However, the corresponding temperature range is reproduced to much bet- 
ter accuracy than that of [20J, where the calculated temperature Trm is far too low. The 
magnetic-susceptibility data follow roughly the Curie- Weiss law with the paramagnetic Curie 
temperature 0m ~ 120 K. Only the dashed curve in Fig. [7] reproduces correctly the approx- 
imate Curie- Weiss law at high temperature. This unusual behavior from a magnetic point 



of view can be understood easily from Eq. (28), where the nonlinearity in magnetization 



(the term oc M 2 ) can be dominated by the magnetoelectric coupling, as Tf/T M » 1. 
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FIG. 6: Dielectric susceptibility as a function of temperature for the specific values of applied 
magnetic field. 



C. Arrot planes 

In systems with one order parameter a convenient way to represent their behavior near 
the phase transition temperature is to draw Arrot plots. For a ferromagnet, the Arrot plot 
is a series of curves M 2 vs. H a /M. The same approach can be used in our case with two 
order parameters. The only difference is that now instead of having a series of curves we 
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FIG. 7: (Color online). Magnetic susceptibility as a function of temperature. (Solid line) The 
inverse molar magnetic susceptibility of BiMnO"3 measured in IT |19] has been rescaled by its value 
at T = 101-KT in order to obtain the dimensionless quantity. Dashed and dotted lines represent the 



respective calculated temperature variations of the magnetic susceptibility ( 28 ) ; we used the same 
parameter values as in Fig. [2} 

have sets of planes for each of the order parameters. Namely, we have the dependencies 
M 2 , P 2 (j^, =^J. The plane which crosses point (0,0) corresponds to the phase 

transition temperature in zero field. A representative set of the Arrot planes is drawn in Fig. 
[8] for the magnetization. One can see that the transition temperature can be determined 
from the dependence M 2 vs e/P, not only from M 2 vs h/M\ 

From (3) and Q after dividing both equations by the order parameter and solving the 
resulting system we obtain: 
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From Eq. ( 29 ) one can also see that for T > Trm we still have a nonzero value of P (cf . 



Eq. (20)). Therefore, the Arrot planes P 2 (H a /M, E a /P) are not important, since Tf > Trm- 
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FIG. 8: (Color online). Arrot set of planes for the magnetization, calculated for BiMnO"3 for selected 
temperatures. The dark plane correspond to the phase transition temperature. The parameters 
are the same as in the earlier Figures. 

D. Specific heat 

In the framework of Ginzburg-Landau theory we may also calculate the specific heat as 
the second derivative of the free energy. 

Providing a similar analysis as in the previous Section we can write down the following 
expressions for the specific heat valid in the respecting temperature regions: 

2(1-73.) b 1 i° T 1 < I (32) 

| • f T f° rT > T RM- J 
Where ao = l/(T^s/jmo), b = 1/(5^). The magnetoelectric part AC P of the specific 
heat calculated in this manner is shown in Fig. [9j The mean field values are close to the 
experimental data below the transition temperature Trm- Essential differences appear above 
the transition and may be attributed to either short-range or fluctuation effects, as discussed 
in the next Section. Nevertheless, in spite of the discrepancy in AC P the Landau approach 
predicts the basic characteristics curve of this magnetoelectric system, surprisingly well at 
low temperatures T < Trm- 




IV. GAUSSIAN FLUCTUATIONS OF THE ORDER PARAMETERS 



The fluctuations of the magnetization seem to have a small effect on the magnetization 
curve close to T RM , but there is a discrepancy for the specific heat AC P (T). AC P (T) is 
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FIG. 9: (Color online). Temperature dependence of the specific heat for BiMnOs: mean-field 
specific heat values - dashed line, experimental data [13] - solid line. The magnetic part of the 
specific heat data are obtained by subtracting the total specific heat of BiSc03 from that of BiMnO"3 
|X3|, I22j . since the former is not magnetic. 

taken as the difference between the total specific heat of BiMn03 and that BiSc03 |13j . In 
this manner, AC P (T) represents only the magnetic part of the specific heat. Therefore, we 
discuss the role of Gaussian fluctuations of M(r) and P(f) on the thermodynamic properties 
next. 

A. Landau functional in spatially inhomogeneous case 

In previous Sections, we made a very crucial assumption namely, that the order parame- 
ters are spatially homogeneous. The following calculations are provided in order to improve 
the temperature dependence of the specific heat part AC P (T) obtained above by taking into 
account the spatial fluctuations of the order parameters. For that purpose, we introduce the 
effective free energy F as a functional of spatially inhomogeneous order parameters [23] : 




where (M, P, T) is the free energy functional which incorporates thermal fluctuation in 
equilibrium, i.e. 





(33) 
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<f>(M(r),P(r),T) = Q ° (T - Tm) M(r) 2 (34) 

+ jM(f) 4 + ^|VM(f)| 2 

+ ^-7>) p(rr + ^ (rl4 

+ ^|VP(f)| 2 + |(P(f)M(f)) 2 . 

In the present situation, we use again dimensionless units and then the functional 
F{M, P, T} takes the form: 

tta 2 n T, 



where: 



Fo + q^M_ Jj(M(r),P(r),T)d\ (35) 



M(v),P(v),t) = |P 2 + |P 4 + | |VP| 2 (36) 



b 11 
+ -M 2 + -M 4 + -|VM| 2 
2 4 2 



We use the following rescaling: r = |, £ = v/^§^> ai = ^ (^w) (T/Tf — 1), a 2 = / 
°3 — n T^"' n = iwh = — 1' c m = Imj^-, and integrate over the volume VL = £ 3 . In 
doing so, we assume that the volume dependence is determined by £ 3 , which will be regarded 
as a fitting parameter. 

B. Effect of Gaussian fluctuations 

The starting point for the following calculations is to consider small fluctuations around 
the mean field values of the order parameters: 

M(r) = M + 8m(r), (37) 

P(r)=P + 5p(r). (38) 
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After substituting (37) and (38) into (37) we expand the expression for the free energy 



density. For simplicity, we retain only the second order terms: 

a l r->2 , °2 ni , b\ 



} (M(r),P(p),r) « |p 2 + f P 4 + |m 2 (39) 

+ \m£ + °fP 2 M 2 

+ 5p ( ai Po + a 2 P* + c m P 2 M ) 

+ 5m (6iM„ + M 3 + c m P M 2 ) 



+ (5m) 2 (| + ^M 2 + i Cm P 2 

+ —P M 5m5p 

+ l -\V5m\ 2 + ^\V5p\ 2 . 

_ _ 



The constant expression: 



0o = |P 2 + | P 4 + |M 2 + ^M 4 + ^P 2 M 2 , (40) 

gives the mean filed value of the free energy density. Linear terms in {5m, 5p} vanish, 
because the mean field solution {Mo, Po} minimizes the free energy F. In effect, we obtain 
the contribution to the free energy coming from the fluctuations of the order parameters in 
the form: 

5F = C m J U m (5p) 2 + B m (5m) 2 + l -\V5m\ 2 + j\V5p\ 2 + lme 5m5p\ d 3 r, (41) 

where: A m = f + |a 2 P 2 + \c m M 2 , B m = | + \M 2 + \c m P 2 , 7me = c -fP M Q , and C m = 
b 

With the help of the expression for the 5F obtained above we calculate in the Appendix 
A the explicit form of the partition function Z, which takes into account spatial fluctuations 
of the two interacting order parameters. 

C. Specific heat 



To calculate explicitly the contribution of the fluctuations to the specific heat we use the 
partition function (A17). The free energy part due to fluctuations can be written as: 
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SF = -k B T\nZ (42) 




P m + ia 3 k 2 )(5 m + ik 2 ) + hi 
We change the summation over k into integration and have: 



A F = JJL [ kmax KfMY 1 \ mo, 

2 A (2tt)3 D UicJ \{A m + §a 3 k 2 )(£ m + |k 2 ) + i 7 LI ? ' 1 



After differentiating twice (43) (cf. Eq. (31)), we obtain the part 5C P for the specific 
heat including Gaussian fluctuations. Finally, the total specific heat becomes 

AC r = -T^-'f + SC p , (44) 

where we have added 5C P to the mean-field part. Following our previous notation we set 
the integration limit as k max = where we take the lattice parameter a = 9.5415A |13j . 



In Fig. 10 we compare the theoretical results for the specific heat with the temperature 
dependence of AC P = C p \BiMn0 3 ~ C p \BiSc0 3 - We see that the fluctuations overestimate 
the experimental behavior for T < T RM and underestimate the data for T > T RM . We 
attribute this (cf. Sec. Ill B) to the role of short-range order which gradually disappears 
as T increases above Trm- This is also the reason why the mean-field results match the 
experiment well for T < Trm- 



D. Correlation lengths in multiferroics 

In an analogous manner one can calculate the evolution of the correlation length through 
the magnetic phase transition. 

The correlation function of an order parameter ip(r) in two distant points is defined as: 

c,(f, P) = ^(V>(r) - - V>)\ = ^5ip(r)5i/;(r')^ . (45) 

After a Fourier transform we can write g(f,r') as: 

(5iP(fi)6iP{f 2 )) = (^fae'^-^e^-*)) (46) 
fei,fc 2 
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ro 80 90 100 110 120 130 

Temperature (K) 

FIG. 10: (Color online). Temperature dependence of the specific heat for BiMnO"3. Dotted line: 
the mean-field part of the specific heat; thick solid line: experimental data |13] was taken as 
C p \BiMn0 3 — Cp\BiSc0 3 - Thin solid line: specific heat after taking into account the Gaussian 
thermal fluctuations calculated for £ = 2.71 • 10 _8 [m] and n = 1. 



and finally: 



(47) 



where r = r~2 — f\ . 

To evaluate the coherence length for both the magnetic and the electric order parameters 



we need to calculate {Sm.^Sm^) and (<5p- k <5p k ). Using (A7) and (A8) we obtain the following 
expressions: 

(5m. k 5m k ) = (5m 2 lk + 5m 2 2k ) , (48) 
Wft) = <^ 2 lk + £j4) • (49) 



(5p. k 5m k + <5p k <5m_ k ) = (5pi k mi k + 8p 2 \jn 2 \ s ) • 



(50) 



The average values (5m. k 5m k ) and (5p. k 5p k ) can be calculated using the matrix A k denned 
in the Appendix A: 



(5m. k 5m k ) 



2 ) 

2k) 



J Vr) V e-^^ A ^l(5p 2 k + 6p 



2k) 



J Pr/ k e~^k'£ A kiC 



(51) 
(52) 
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If Ak is symmetric for real Gaussian integrals as in (51) and (52) we have: 



where Ay - is an element of the Ak matrix in the zth - row and jth-column, rji for % 
is an element of the vector ffk. 
We obtain following expressions: 



(<WPk> = Ail + A 



1 A-l _ k B T 



2C rn {A m + \a z k 1 ) 
k R T 



(5m. k 5m k ) = A 3 l + A J = - , .,, ■ 



2 

and for the cross-correlations we have: 



(6pi k rai k + 5p 2 k"i 2k ) = A 3 l + Aj 



k B T 



Cm'Jme 

Finally, the correlation functions for the magnetic and electric subsystems take 



g P {f[,r 2 ) = (5p{fi)5p{f 2 )) 

d 3 k 



(5p. k 5p\l} e 



ik-(r"2— rl) 



(2tt) 3 

,ik-(r"2— fx) 



d 3 k k B T 
~Wf ' 2C m (A m + U 3 k 2 " 



gm{fuf 2 ) = (5m(f[)5m(r 2 )) 



d 3 ^ _ k B T ik-(r2-rl) 



Using: 



(2vr)3 2C m {B m + \k 2 ) 
d 3 k 1 



ik-f 



(27r) 3 k 2 + a 2 47rr 
we obtain the correlation function in the Ornstein-Zernike form: 



^ v fc B T e V -3 |r2 ri1 fc B T e 

g P {ri,r 2 ) = — — — = — — — 

C m a 3 47r|r 2 -ri| C rn a 3 Aix\r 2 - r x \ 
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k R T e -V2B^|r5-fi| kl}T 



gm[ri,r 2 ) 



with the correlation lengths: 



e 



C m 47r|r2-rl| 



Q3 



1 



2S„ 



(60) 



(61) 



From Eqs. (59) and (60) we see that when the coherence length is close to zero the 
correlation function becomes equal to zero as well. On the other hand when the coherence 
length is large the correlation function decreases as oc , ^ 1 The correlation radii in this 
region are significantly greater than the lattice constant. 

2r 



— i 

bC 

U 
fH 

U 



0.5 





■ ■■■ /= 



















200 400 600 
Temperature (K 



500 



FIG. 11: (Color online). Temperature dependence of the coherence lengths for magnetic (solid 
line) and electric (dotted line) subsystems calculated for BiMnO"3. The parameters we use are the 
same as in the earlier Figures. 



We can see that close to ferroelectromagnetic phase transition temperature both coher- 
ence lengths behave in a similar way. Hence, we assume that below Trm, where both 
magnetic and electric ordering is observed, the coherence length coalesces as T decreases 



(T — > 0) (cf. Fig. 11). It can also be seen (from the Fig. 11) that the order parameter cor- 
responding to the lower critical temperature experience smaller fluctuations near the phase 
transition. 
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We calculate the cross correlations of the order parameters, which take the form: 



9pm(fi,r 2 ) = (Sm(fl)Sp(r 2 )) (62) 
= ^2 ($ m ik5pik + 5m 2 kSp2k) e lk r 

k 

1 k B T 



2tt 2 C m ^ me r Jq 
k B T 



k sin(kr)dk 



. , , 3 (svd-(k max r) T'k max cos(k max r)) . 

We see that the correlation length £ p near and above the magnetic phase transition is 
enhanced, whereas the cross correlations exhibit an oscillatory behavior and evolve contin- 
uously through Trm- 



V. CONCLUSIONS 



In this paper we examined a simple Landau approach for a system with two single- 
component order parameters representing ferroelectricity and ferromagnetism, respectively. 
The obtained results are consistent with experiment. Magnetoelectric coupling introduces a 
strong renormalization of the ferromagnetic transition temperature (enhancing it by a factor 
of 50) causing the magnetic phase transition to be observed at T RM = 100.5 K. A number of 
coupling effects such as the enhancement of electric polarization, the anomaly in the dielectric 
susceptibility, and a fairly large negative magnetocapacitance, occur concomitantly in the 
vicinity of T RM . We introduced a simple extension of the Arrot plot which we called in 
text Arrot planes. We consider this concept to be potential useful in evaluation of the 
phase transition temperature while using only field dependence of unbounded with it order 
parameter. We have also estimated the contribution of the Gaussian fluctuations to the 
specific heat and have noted that short-range-order effects are not accounted properly. It 
would be important to understand the microscopic reasons of such a strong enhancement of 
the magnetic exchange interactions caused by a monoclinic lattice distortion which leads to 
the appearance of the ferroelectric dipole moments. 
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Appendix A: Partition function with inclusion of Gaussian fluctuations 

If we assume that the fluctuations do not change rapidly in space then we can estimate 
their local value in terms of their Fourier components: 



k 

5m(r) = <5m k e 



ikr 



where k = £k. We also assume that e lk ' r satisfies periodic boundary conditions. 



After substituting (Al) and (A2) into (41): 



(Al) 
(A2) 



k q 



<5m k 5m q / d 3 re i(k+q) r 



5F = C mlme J2J2 5m ^ I d 3 re i( - k+ «> r (A3) 

k q 

+ C m ^ ] ^ ] B rn — — k ■ 

5p*5 P(l J d 3 re* (k+q)r . 

Taking into consideration that f rf 3 re^ k+q - ) r = <5 krq and after rewriting the last term as 



k q 



A m - 2 a 3 k ■ q 



k q 



- Y [^ m k^P-k + Sm.^Spk] , 



(A4) 
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we obtain the following expression for the free energy: 

5F = C m J2( A ™ + ^a 3 k 2 )5p k 5p. k (A5) 

k 

+ C m ^2(B m + ^k 2 )5m k 5m. k 

k 

+ T^Cmlme ^ [^ m k^-k + Sm. k 5p k ] . 
k 

In order to evaluate the value of the fluctuations we have to take an average of all the 
possible configurations. One can define the statistical sum for a system with two order 
parameters as an integral over all existing profiles for each of the order parameters: 

Z = JJ / V (5m k ) V (Spk) e *b* . (A6) 
k J 

Because 5m k and Sp k are complex numbers we can represent them in the following way: 

5m k = 5m lk + i ■ 5m 2k , (A7) 

5p k = 5p lk + i ■ 5p 2k . (A8) 

As may be seen 5m k and Sm. k are not independent. In order to avoid double counting for 
5m k and 5m. k we have to take only wave vectors k with k z > (the same applies for 5p). 
Thus, the partition function representing the Gaussian fluctuations takes the form: 

r°° r°° r°° r°° sF(s m ,s P ) 

Z= j [ / V(5m lk ) V(5m 2k ) V (Sp lk ) V (5p 2k ) e ^ , ( A 9) 

kfc >0 00 J —oo J — oo 

where the corresponding free energy functional is 

5F = 2C m ( A ™ + l^ 2 )(5p 2 lk + 5p 2 2k ) (A10) 

k,k z >0 

+ 2C m {B m + l -\t 2 ){5m 2 ^ + 5m 2 ^) 



k,fc z >0 

+ 2C m 7 me ^ (^lk^pik + ^makfek)- 

k,fc z >0 



For further calculations we use matrix notation: 
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Vk = 

\5m 2 kJ 
tfl = {Spik, Sp 2 k, Sm lk , 5m 2k ) 



(All) 



(A12) 



and 



+ \ci3k 2 



Al 





2Cm7n 







AC m (A m + |a 3 k 2 




2C m 7 me \ 

2,C m 'y me 
AC m (B m + ±k 2 ) 

AC m (B m + Ik 2 ) J 



(A13) 



In these terms the statistical sum Z can be written as: 



J l. 



(A14) 



Because the matrix Ak is symmetric and the vector rfk is real, we can use the expression for 
the last Gaussian integral: 

Z k = (27r)£(de*A k )-3, (A15) 
where D is the dimension of the vector fj k which in our case is equal to 4. The final expression 



for Zk after diagonalization of (|A11|) can be written as: 

1 



TV 2 (k B T 



' k 4 V C m J \{A m + \a,k 2 ){B m + \k 2 ) - i 7 ^ e| 
The corresponding total statistical sum is: 



(A16) 



n * 



z = i i ^ k 

k,fc z >0 



(A17) 



This expression is used in Sec. IVC to calculate the specific heat. 



Appendix B: Angular degrees of freedom in order parameter fluctuations 
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In our calculations we see a big discrepancy in specific heat 
near the phase transition Tm- Here we present a brief discussion 
of possible way to improve accordance to the experimental data 
by considering a little bit more realistic model. We assume that 
our order parameters are three dimensional. Due to the fact that 
our main interest lies in region near ferromagnetic phase transition 
and because Tf » T M we can consider a following simplification: 
FIG. 12: (Color online). P(f) = P{r)e z . The main idea of proposed improvement is to take 
Model picture of the or- ' lT iio account the angular fluctuation of M near the phase transition, 
der parameters arrange- If we apply the orientation of the order parameters as presented at 




ment. Fig. 12 we obtain following Landau free energy potential 



AF = a -M 2 z cos 2 6 + h -Mt cos 4 9 + |p 2 + ^P 4 + 7 (M 2 ) 2 cos 2 6. (Bl) 

In the spirit of mode-mode coupling approach we can present the following two terms 
from IBll as: 



(M z cos(9)) 2 « (M 2 ) cos 2 6 + (cos 2 6) M 2 - (M 2 ) (cos 2 0) (B2) 

and 

(M z cos(0)) 4 « {M 2 Z ) (cos 2 6) {Ml + (M 2 ) cos 2 9 - (cos 2 6) (M 2 )} . (B3) 



After substitution B2 and B3 in Bl the free energy potential will have the form: 



AF = Ml (cos 2 6) + \ (Ml) (cos 2 6) + |p 2 (cos 2 6)} (B4) 
+ cos 2 6 { I (Ml) + \ (Ml) 2 (cos 2 6) + |P 2 (Ml) 
+ P 2 {| + ^(P 2 > + |<M 2 ><cos 2 e>} 



(M 2 2 ><cos 2 e>{V + ^(M 2 ><cos 2 e> 
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